Preliminary results for RR Lyrae stars and Classical 
Cepheids from the Vista Magellanic Cloud (VMC) Survey 
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Abstract The Vista Magellanic Cloud (VMC, PI M.R. 
Cioni) survey is collecting Ks-hand time series photom- 
etry of the system formed by the two Magellanic Clouds 
(MC) and the "bridge" that connects them. These data 
are used to build i^s-band light curves of the MC RR 
Lyrae stars and Classical Cepheids and determine ab- 
solute distances and the 3D geometry of the whole sys- 
tem using the A'-band period luminosity (PLKs), the 
period - luminosity - color (PLC) and the Wesenhiet 
relations applicable to these types of variables. As an 
example of the survey potential we present results from 
the VMC observations of two fields centered respec- 
tively on the South Ecliptic Pole and the 30 Doradus 
star forming region of the Large Magellanic Cloud. The 
VMC i^5-band light curves of the RR Lyrae stars in 
these two regions have very good photometric quality 
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with typical errors for the individual data points in the 
range of ^ 0.02 to 0.05 mag. The Cepheids have ex- 
cellent light curves (typical errors of ~ 0.01 mag). The 
average Ks magnitudes derived for both types of vari- 
ables were used to derive PLKs relations that are in 
general good agreement within the errors with the lit- 
erature data, and show a smaller scatter than previous 
studies. 

Keywords stars: variables: Cepheids- stars: vari- 
ables: RR Lyrae ~ galaxies: individual: LMC - galax- 
ies: distances and redshifts 



1 Introduction 

The VISTA near- infrar ed YJKc; survey of the Mag- 
ellanic system (VMC; ICioni et all I2OI1I ) is an ESQ 
public survey that is obtaining deep near infrared 
imaging in the F, J and Ks filters of a wide area 
acro ss the Magellanic sy stem, using the VIRCAM cam- 
era dOalton et all 120061) of the ESQ VISTA telescope 



( Emerson. McPherson. fc Sutherland 1120061 ) 



The main science goals of the survey are the deter- 
mination of the spatially-resolved star-formation his- 
tory (SFH) and the definition of the three-dimensional 
(3D) structure of the whole Magellanic system. The 
VMC observations are devised to reach Ks ^ 20.3 mag, 
thus allowing us to measure sources encompassing most 
phases of stellar evolution: from the main-sequence, to 
subgiants. upper and lower red giant branch (RGB) 
stars, red clump stars, RR Lyrae and Cepheid vari- 
ables, asymptotic giant branch (AGB) stars, post-AGB 
stars, planetary nebulae (PNe), supernova remnants 
(SNRs), etc. These different stellar populations will 
help us assessing the evolution of age and metallicity 
within the whole MC system. 
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111 this context, a very significant role is played by 
the radially pulsating variables, since these stars obey 
a period-mean density relation that forms the basis of 
their use as standard candles to measure the distance 
to the host system. The RR Lyrae stars in particu- 
lar, in the K band obey a period-luminosity-metallicity 
(PLZ) relation which is weakly affected by the evo- 
lutionary effects, the spread in stellar mass within 



ing corrections (see e 


.e. Lonemore. Fernlev. & Jameson 


19861 


DairOra et al. 


2004; CoDDola et al. 


2011). Sim- 


ilarly, the Cepheid PL relation in the 


K band is 



much narrower than the corresponding optical re- 
lations, and less affected by systematic uncertain- 
ties in the reddening and meta l content (see e.g. 
Caputo. Marconi, fc MusellalboOOt) . 



We present preliminary results for the RR Lyrae 
stars and Cepheids contained in the first two "tiles" 
completely observed by the VMC survey, namely tiles 
8_8 and 6_6. These two tiles are centered respectively 
on the South Ecliptic Pole (SEP) field and on the well 
known 30 Doradus (30 Dor) star forming region of the 
Large Magellanic Cloud (LMC). The SEP field is par- 
ticularly interesting because this region of the sky will 
be continuously and repeatedly observed during the the 
commissioning phase of the astrometric Gaia satellite 
just after the launch in Spring 2013. 



2 The VMC data for the variable stars 

The VMC o bserv ing strategy is described in detail in 
Cioni et al.l (j201l[ ) to which the interested reader is re- 
ferred for more information. Here we focus on the data 
acquisition for the variable stars. 

In order to obtain well sampled light curves and mea- 
sure accurate parameters (namely Ks average magni- 
tudes) for the variable sources, the Ks observations 
were split into 12 epochs and distributed along several 
consecutive months. 

We used the vl.O VMC release "pawpriiits" (6 "paw- 
prints" form a "ti le"). The pawprint s were processed 



by the pipeline (llrwin et al. I 120041) of the V ISTA 



Data Flow System (VDFS. lEmerson et al. \\2004 and 
retrieved from the VISTA Science Archive (VSA, 
Hamblv et aDl2008ll^ . 



Usually a variable star is present in two or three not 
necessarily consecutive pawprints. Hence we first cal- 
culated a weighted average of the pawprints' Ks mag- 
nitudes to obtain the "tile" Ks, which also represents 
one epoch of data. 
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Fig. 1 A's-band light curves of three fundamental mode 
(F; 3 upper panels) and one first overtone (FO; bottom 
panel) RR Lyrae stars in the SEP field. Internal photo- 
metric errors of the individual data points are shown, they 
are typically in the range of 0.02 to 0.05 mag. Periods were 
taken form the EROS-2 Survey. 



■"-http:/ /horus. roc. ac.uk/vsa/ 
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Fig. 2 Ks-h&nd light curve for a sample of Cepheids in the 
SEP field. Note in the third panel the nice light curve for 
a faint FO Cepheid. Periods were taken form the EROS-2 
Survey. 



Particular care was devoted to the determination of a 
proper Heliocentric Julian Day (HJD) for the Ks value 
of each "tile" . 

The identification, pulsation period and epoch of 
maximum light of the RR Lyrae stars and Cepheids 
contained in the SEP and 30 Dor regions are already 
available in the photometric archives of the Microlens- 
ing surveys OGLE and EROS-2. 

The OGLE surveys (jSoszvnski et al.lboosllioogl and 
references therein) of which stage IV is in progress, 
cover an area of the MC system that extends pro- 
gressively further outside from the bar of each of the 
Clouds. EROS-2 (|Tisserand et al.lboOTi) is, at present, 
the most extended of the two surveys and is the only 
one covering the largest fraction of the VMC field of 
view including the most peripheral areas such as the 
SEP region, that in fact only overlaps with EROS-2. 
Coordinates, periods and optica^ light-curves of the 
RR Lyrae and Cepheids in the Gaia SEP field were 
taken from the EROS-2 catalogue and cross-matched 
to the VMC catalogues for the SEP tile. The 30 Dor 
field is covered by both EROS-2 and OGLE-III, but for 
the present analysis we only employed the periods and 
optical (Jo hnson-Cousins V, I bands) l ight curves from 
OGLE-III ISoszvnski et al.ll2008i liooj ). 
As a result of the matching procedure we found 107 
and 1465 RR Lyrae variables, and 11 and 326 Classical 
Cepheids in the SEP and 30 Dor fields, respectively. 

The periods available from EROS-2 for the SEP vari- 
ables and from OGLE-III for the 30 Dor variables were 
used to fold the ifg-band light curves produced by the 
VMC observations and derive average Ks magnitudes 
for the variables. Examples of the VMC Ks-hand light 
curves of RR Lyrae stars and Cepheids in the SEP and 
30 Dor regions are shown in Figs. [U [5] and [31 |31 re- 
spectively. Error bars of the individual Ks measure- 
ments are shown in the figures. The uncertainties for 
the Cepheids are of the same size at the data-points, 
while they are larger for the RR Lyrae stars, particu- 
larly in the 30 Dor region, due to the fainter magnitudes 
of these variables and the high crowding of the 30 Dor 
field. 



2.1 Determination of the average Ks magnitudes 

In the present work we confine our discussion to the RR 
Lyrae stars in the SEP field and the Cepheids in the 30 
Dor region only. This is because of the small number 
of Cepheids in the SEP tile (only 11), and because the 



^The EROS-2 blue channel (420-720 nm), overlaps with the 
V and R standard bands, and the red channel (620-920 nm) 
roughly matches the mean wavelength of the Cousins / band 
jTisserand et al.ll200'ii) 
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high photometric contamination at the magnitude level 
of the RR Lyrae stars in the 30 Dor field requires a more 
careful and detailed analysis which is still in progress. 

The A'5-band light curves of the RR Lyrae stars 
are almost sinusoidal and the amplitude of the lumi- 
nosity variation in Ks is also rather modest. How- 
ever, the most accurate results for the mean magni- 
tudes are ac hieved by u s ing th e light-curve templates 
provided by iJones et al. I d 19961) . Note that the use of 
templates requires accurate ephemeris and l^-band am- 
plitudes which, for the SEP RR Lyrae were calculated 
starting from the EROS-2 database. 

As for the Cepheids, as shown by Figs. [21 21 the 
light curves are well sampled and nicely shaped. Hence, 
to derive the (intensity) average < Ks > we simply 
used a custom program written in "c" language that 
performed a spline interpolation to the data. 



3 Results 

Having derived the mean Ks magnitudes as discussed 
above, we then computed the PLKs relations for the 
SEP RR Lyrae stars and the 30 Dor Classical Cepheids. 

3.1 RR Lyrae stars in the SEP field 

The PLKs relation of the 95 RR Lyrae stars with good 
photometry in the SEP field is shown in Fig. [5l where 
the periods of the FO RR Lyrae stars were fundamen- 
talized using the formula logPi?=logPFO + 0.127. The 
Ks magnitudes were dereddcned using the E{B — V) 
value for the SEP field provided by S. Rubele (pri- 
vate communicatio n) and calculated as described in 
Rubele et all |201l[ ). 

We then performed a least-square fit to the data, ob- 
taining the following PLKs relation: 



Kl 



(17.351 ± 0.035) - (2.450 ± 0.133)Zo5P 



(1) 



which has an r.m.s. =0.082 mag. The higher dis- 
persion of the SEP RR Lyrae PLKs with respect 
to the PLKs relations followed by the RR Lyrae in 
Galactic globu l ar clu sters (see, e.g. DaH'Ora et al.lliooi : 
Coppola et al.l I2OIII and references therein) is likely 
due to a spread in distance among the SEP RR Lyrae 
along the line of sight, as well as to the possible pres- 
ence of a metallicity variance among the variables. We 
have compared our results for the SEP RR Lyrae stars 
with previous near-infrare d observations o f RR L yrae 
variables in the LMC by iBorissova et al.l (j2009[ ) and 
Szewczyk et al.l (j2008[ ). Our Eq. [T] can be compared 



with Eq. 7 of IBorissova et al.l (|2009f ) which was calcu- 
lated using RR Lyrae data from both the above men- 
tioned papers (107 stars in total): 



K% 



(17.47 ± 0.07) - (2.11 ± Q.n)logP 



(2) 



The two equations a re only marginally co nsistent 
within the errors, and the Borissova et al. ( 20091 ) PLKs 
has a larger dispersion than ours. These occurrences 
can be due to both geometrical (depth) effects and a 
metallicity spread. These aspects will be investigated 
in a following paper (Moretti et al. in preparation) 
where we shall also provide our estimate of the LMC 
distance. 

3.2 Classical Cepheids in the 30 Dor Field 

The PLK relations of the 172 fundamental mode and 
the 154 first overton^ Classical Cepheids analysed in 
the 30 Dor field are shown in Fig. |6l Since the satura- 
tion level of the VMC survey in A'5 limits the length 
of the periods we are able to measure to a maximum 
value of about 15-20 days, we decided to c omple ment 
our data with the sample by Persson et al. ( 20041 ) that 
includes 84 F-mode Cepheids with periods mainly rang- 
ing between 10 and 100 days. To merge the two samples 
we first transformed Persson et al.'s original photome- 
try fror n the LCO to the 2MASS system using the rela- 
tions bv lCarpenterl (j200ll ) These data are shown as light 
blue circles in Fig. |6l To account for the variable red- 
dening that characterizes the 30 Dor field, we adopted 



the re cent evaluatio ns bv [ Haschke. Grebel. fc Duffau 



( 20111 ). while for the iPersson et al.l (|2004l ) data set we 
adopted the values provided by the Authors. Finally, 
wc performed least-square fits to the data of F- and 
FO-modcs variables separately, adopting an equation 
of the form K'g=a+P logP. The coefficients derived 
from the fits are provided in Table [H The same ta- 
ble also shows a coinpariso n with the literature re- 
sults of iPersson et al.l (|2004l ). who only published the 
PLKs relation of fu ndamental-mode pulsators, and of 
Groeneweeen I (I2OOOI ) who also gives the FO PLKs re- 
lation. An inspection of the table reveals that 1) for 
the F-mode pulsators there is general agreement within 
the errors among the three studies; 2) the present rela- 
tion has smaller errors of the coefficients of the PLKs 
equation because of the wide range of periods spanned 
by the Cepheids used to derive the relation, including 
for the first time in the near-infrared range a significant 
number of objects with periods smaller than 5 days; 3) 



•^Double mode pulsators F/FO and FO/SO were included in the 
F and FO samples, respectively. 
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Table 1 Jfs-band PL relations for F and FO Classical Cepheids. The model equation is: Kg^a+P logP. The first two 
lines of the table show the coefficients of the least square fit to the data of F and FO Classical Cepheids in the 30 Dor 
region presented in this paper. The last three lines show the coefficients of the same relations found in the literature. 



mode 


Q 




/3 


op 


r.m.s. 


source 


F 


16.070 


0.017 


-3.295 


0.018 


0.102 


This work 


FO 


15.580 


0.012 


-3.471 


0.035 


0.099 


This work 


F 


16.051 


0.05 


-3.281 


0.040 


0.108 


Persson et al. f 20041 


F 


16.032 


0.025 


-3.246 


0.036 


0.168 


Groenewesen (2000) 


FO 


15.533 


0.032 


-3.381 


0.076 


0.137 


Groenewegen (2000) 



17.50 



_ p= 0.4736372 d Mode=F 

* * I i i * ' * I i i * I* 



P= 0.5702073 d Mode = F 



i i I 



iii 



_ P= 0.5916191 d Mode = F 



17.50 - J i 



18.00 



P= 0.5610799 d Mode=F 



18.00 



18.50 



1.0 
PHASE 



Fig. 3 Same as Fig. [T] for RR Lyrae stars in the 30 Dor 
region. Periods were taken form the OGLE-III Survey. 
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Fig. 4 As in Fig [2] but for Cepheids in the 30 Dor field. 
Typical errors of the individual data points are of the order 
of 0.01-0.02 mag, hence have same size as the data points 
plotted in the figure. Periods were taken from the OGLE-III 
Survey. 
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for the FO pulsators the agreement with iGroenewegen 
is not as good as for the F-mode Cepheids. This 
could be due to t he fact tha t the deeper photometry of 
the VMC survey |Groenewegen ,200d rehes on 2MASS 
and DENIS data) allowed us to reach the fainter FO- 
Cepheids that populate the low-period tail of the PLKs 
relation. 

The use of the PLKs relation derived in this pa- 
per to estimate the distance of the LMC as well as the 
discussion of the Wesenheit and PLC relations is de- 
ferred to a forthcoming paper (Ripepi et al. 2012, in 
preparation). 



4 Conclusions and future perspectives 



-0.3 
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Fig. 5 (ifs)-band Period-Luminosity relation for RR 
Lyrae variables in the SEP field. Open and filled circles 
show first overtone and fundamental mode RR Lyrae stars, 
respectively. The solid line shows the result of the least 
square fit to the data (see text for details). 



30 Dor + Persson et al. (2004) sample 




We have presented preliminary results on the RR Lyrae 
stars and the Classical Cepheids observed by the Vista 
Magellanic Cloud (VMC) survey in two fields of the 
LMC, centered on the South Ecliptic Pole (SEP) and 30 
Dor region, respectively. The Ks light curves of the RR 
Lyrae variables have very good photometric quality, for 
the brighter Cepheids we obtained excellent data. The 
average Ks magnitudes derived for both kind of ob- 
jects allowed us to derive PLKs relations which were 
compared with the literature, finding in general good 
agreement within the errors. Our data show in any 
case a smaller scatter with respect to previous studies. 
We expect much better results when more data from 
the VMC survey will be available. In particular we will 
be able to use the PLKs relation for the RR Lyrae vari- 
ables and the PLKs-, PLC and Wesenheit relations for 
Classical Cepheids to investigate both the 3D geome- 
try and the absolute distances to the Magellanic Cloud 
system. 
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Fig. 6 Ks-haxiA PL relation for Cepheids in the 30 Dor 
field. Black open and filled circles show FO- and F-mode 
pulsators, respectively. Li ght blue filled circles show the F- 
mode Cepheid sample bv lPersson et al] l)2004l ). Solid lines 
are the result of the least square fit to the data (see text for 
details). 
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